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ABSTRACT 

We study the one dimensional partially asymmetric simple exclusion pro- 
cess (ASEP) with open boundaries, that describes a system of hard-core par- 
ticles hopping stochastically on a chain coupled to reservoirs at both ends. 
Derrida, Evans, Hakim and Pasquier [J. Phys. A 26, 1493 (1993)] have shown 
that the stationary probability distribution of this model can be represented 
as a trace on a quadratic algebra, closely related to the deformed oscillator- 
algebra. We construct all finite dimensional irreducible representations of 
this algebra. This enables us to compute the stationary bulk density as well 
as all correlation lengths for the ASEP on a set of special curves of the phase 
diagram. 
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1 Introduction 



The one dimensional asymmetric simple exclusion process (ASEP) has been 
extensively studied by mathematicians and physicists (see p| and refer- 
ences therein) as one of the simplest model of a system out of equilibrium. It 
is a diffusion model with hard-core exclusion that can be used to describe hop- 
ping conductivity in superionic conductors, traffic flow, and interface growth; 
it can be mapped to a lattice version of the Kardar-Parisi- Zhang equation 

Subject to open boundaries, where particles are injected or extracted, 
the ASEP displays a rather rich behaviour; it exhibits phase transitions in 
the thermodynamic limit. The exact steady state for the fully asymmetric 
case (particles jump only in one direction) was computed in M for special 
choices of the input parameters and more generally in [H ana^ These 
computations rely on a recursion in the system size. An elegant way to 
exploit the recursive property of the steady state is the matrix Ansatz used 
in [^1, where the stationary probabilities are expressed as matrix elements of 
products of operators, which represent particles and holes. 

The same technique was used in for the partially asymmetric exclusion 
process with open boundaries. The operators used here generate a quadratic 
algebra which is simply related to the so-called deformed oscillator-algebra. 
The phase diagram was derived in the thermodynamic limit, and the steady 
state current in each phase was computed. However, since suitable represen- 
tations of the quadratic algebra were lacking general equal-time correlation 
functions in the stationary state could not be calculated. In a recent pa- 
per, Essler and Rittenberg studied Fock representations of the general 
quadratic algebra (that can generically be mapped on the deformed oscilla- 
tor algebra). They first found infinite dimensional matrices, and then gave 
conditions under which these infinite matrices have an invariant subspace of 
finite dimension. These constraints were explicitely written in terms of the 
parameters of the partially ASEP in the case of a one dimensional and a 
two dimensional invariant subspace. One and two dimensional representa- 
tions were constructed, and correlation functions determined with the help 
of these representations. 

In this paper we use standard methods of linear algebra to classify and 
construct all irreducible finite dimensional representations of the deformed 
harmonic oscillator algebra, without having to start with infinite dimensional 
matrices. We show that there is exactly one irreducible representation for 
any finite dimension n, and we find explicit constraints on the parameters 
of the ASEP that make this representation compatible with the boundary 
conditions (cf. Eq. (PP])). These constraints define the locus of points in the 
phase diagram of the exclusion model that are accessible by finite dimensional 
representations. In these regions of the phase diagram we compute the bulk 
density in thermodynamic limit and all the correlation lengths. Our results 
prove some conjectures raised in P|. Morever, to the best of our knowledge 
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the finite dimensional representation presented here have not been reported 
in the hterature.Q The matrices obtained can also be useful in studying 
algebras associated with more general reaction-diffusion processes. 

This paper is organized as follows. In section 2, we recall the matrix tech- 
nique, relate it to the deformed harmonic oscillator algebra, and construct 
all irreducible finite dimensional representations of this algebra. Section 3 is 
devoted to the calculation of the bulk density and the correlation lengths in 
the thermodynamic limit in the regions of the phase space that are accessible 
by finite dimensional representations. The last section concludes with some 
remarks on possible generalizations. Some general properties of the quadratic 
algebra which are shared by all representations are proved in appendix A; 
these properties are used in appendix B to prove the non-triviality of the 
representations obtained in section 2. 



2 Finite dimensional representations of the 
quadratic algebra 

2.1 Matrix Ansatz for the ASEP with open boundaries 

We consider the one dimensional partially asymmetric exclusion process with 
open boundary conditions. Each site i {1 < i < L) of a. one dimensional 
lattice of L sites is either occupied by a particle (xj = 1) or empty (tj = 0). 
The system evolves according to a stochastic dynamical rule: during each 
infinitesimal time step dt, the transitions allowed for the bond {i, i + 1) with 
1 < i < L — 1 are 

10 01 with rate 1 

01 10 with rate x. (1) 

The parameter x is positive and measures the strength of the driving field; 
one can assume with no restriction that x < 1. The model studied here is 
not totally asymmetric, therefore a; > 0. Particles are injected at sites 1 and 
L with rates (1 — x)a and (1 — x)6 and extracted with rates (1 — x)7 and 
(1 — x)P respectively, where 7 and 6 are strictly positive numbers. 
It was shown in |§] that the quantities /(ri, . . . , tl) defined as 

/(ri,...,r^) = {W\f[{T,D + {l-n)E)\V) (2) 

i=l 

^ A reason for this fact may be that, as shown in ||^, there are no irreducible repre- 
sentations of finite dimension, strictly greater than 1, of the algebra associated with the 
totally asymmetric exclusion process, the model mostly studied. 
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provide stationary solutions of the master equation of the ASEP if the oper- 
ators Z), E satisfy the algebra 

DE-xED = {l-x){D + E), (3) 

and the vectors {W\ and \V) are such that 

iPD - 5E)\V) = \V) , 
{W\{aE - jD) = {W\ . (4) 

The ASEP on a finite lattice has a unique stationary state (Perron-Frobenius 
theorem, see e.g. |^). Therefore, all the quantities /(ti, . . . ,tl) are propor- 
tional to the steady state probabilities of the system. If the matrix elements 
(0) are not all equal to zero, then the normalized stationary probabilities of 
tne exclusion process are given by: 

P{ti,...,tl) = -^fin,. . . ,Tl) , 

with Zl = {W\{D + E)^\V) . (5) 

To compute the stationary probability distribution of a system with L sites, 
one needs to find representations of the algebra (H) and of the boundary 
vectors (^). These representations must be such that the matrix elements of 
length L, i.e. of the type 

(VF|L)"^E'"i . . . D''^E""'\V) with + rui + . . . + Uk + ruk = L , (6) 

do not identically vanish (i.e are non-trivial). In appendix A, we find general 
properties, that a representation must possess to ensure that matrix elements 
of the type (0) are non-trivial for a system of length L. 

In [If, relations (|^) and (|) are interpreted as a Fock representation of 
the quadratic algebra generated by the operators 

A = pD-6E -1 and B = aE--fD-l. (7) 

The matrices A and B act trivially on the boundary vectors: 

A\V) = and {W\B = 0. (8) 

Thus \V) and {W\ play the role of right and left vacuum for the operators 
A and B. If a3 =^ 70 the operators D and E can be expressed as linear 
combinations of A and B, which hence satisfy the quadratic algebra: 

(1 - x)a-rA'^ + {a(3 - x^6)AB + (7^ - xa(3)BA + (1 - x)6(3B^ = 
(1 - x) {{a(3 - 75)(a + 7) - a{f3 + l) - 7(« + ^)) ^ 
+ (1 - x) ((a/3 - -fS)iP + S)- SiP + 7) - p{a + 5)) B 
+ (l-x)((a/3-7<5)(a + /3 + 7 + 5)-(a + 5)(/5 + 7)) . (9) 

Conversely, a Fock representation of a general quadratic algebra can gener- 
ically be transformed to a representation of an algebra of the type (^j with 
boundary conditions like (H). 
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2.2 Classification of irreducible representations 

We show that the n dimensional irreducible representation of the algebra (|^) 
can be written in a suitable basis as: 



D 



and E 



( 1 + a 





1 




V 





1 + ax 






1 + ^ 

ax 

1 






1 + ax^ 






1 + - 















1 + ax"-^ / 



1 ^"r-^J 



(10) 



(11) 



a being a non-zero real parameter. 

In order to derive this result, we start by defining two operators d and e 
such that 

D = l+ rf and E = l + e. (12) 

From Eq. (|3|), we see that d and e satisfy the deformed harmonic oscillator 
algebra [p^, |11||: 

de — xed = 1 — x . (13) 
Using equation (|13D, we note that if A is an eigenvalue of the operator de, 

then I 1 H ) is an eigenvalue of the operator ed (we recall that x 7^ 0). 

V ^ J 

In finite dimension the matrices de and ed have the same spectrum, which 
contains therefore the numbers: 



A,l + 



A 



X 



x^ 



1 + 



A- 1 



keN. 



x" 



But there is only a finite number of distinct eigenvalues in a finite-dimensional 
representation: this implies that either x is a root of unity or A = 1. Only the 
latter case is possible here: the eigenvalues of de must all be equal to 1. The 
operator de is invertible, and so are the matrices d and e (their determinant 
can not be zero). We rewrite now equation (13) as: 



d{e — d ^) = x{e — d ^)d . 



(14) 
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Hence, our problem is reduced to finding representations of tlie algebra 

VE = x£T> with V invertible, 

where V = d and £ = {e — d^^) . 

If I a) is an eigenvector of V with eigenvalue a (which is different from 
because D is invertible), then S\a) is an eigenvector of V with eigenvalue 
xa or is the null vector. Consequently, the space spanned by the linearly 
independent vectors {\a),S\a), . . . ,£^\a), . . .} is stable under both V and S. 
To obtain a finite dimensional representation, there must be an integer n such 
that S^\a) = 0. This integer n is also the dimension of the full representation 
space, since the representation is irreducible. 

The proof that the irreducible n dimensional representation of the quadra- 
tic oscillator algebra is given by ([T0|) and (pUf ), is completed by writing 
the matrices oiD = l + d= l + V and E = l + e= l + S + V'^ in 
the basis (la), £|a), . . . , £^"~^|a)). (Hereafter, this basis will be denoted by 
(|l),...,|n)).) 

2.3 Boundary vectors and conditions on the rates 

We construct the boundary vectors \V) and {W\ associated with the irre- 
ducible representation of dimension n found above. These boundary vectors 
are right and left eigenvectors of the operators A and B, defined by Eqs. (|^), 
with zero eigenvalue (see Eq. (||)). In the basis (|1), . . . , |n)) the matrices A 
and B are bidiagonal, their eigenvalues are readily obtained. One can check 
that for A and B to have as an eigenvalue, there must be two integers k 
and I between and n — 1 such that: 

p(l + ax^)-6(l + ^)-l = 
ax'' 

and a(l + ^)-7(l + ax')-l = 0. (15) 
ax' 

Necessarily | / — A;] = n — 1, otherwise a representation of dimension less than 
n would suffice. In the case k = n — 1 and / = 0, one has 

\V) = \n) and {W\ = (1| . 

With such a choice, one has = for any L. This implies, using 

the definiteness property (proved in appendix A), that all the quantities 
/(ti, . . . , tl) are equal to 0: this case has to be excluded. 
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The only case that remains is A; = and / = n — 1; the two equations 
corresponding to (jl^) are 



a 



pa'^ + {/3-5-l)a-5 = 



and they must have a common root a. Both equations have the same struc- 
ture, their solutions can be written a = n±{P,S) and — ^ = K.±{a,'j) with 



^±{u,v) = ^ ^ . (17) 

This function already appeared in the study of the phase diagram of the 
partially asymmetric exclusion process in |^. One can check that for u and 
V positive, 

/t+(-u,f)>0 and -1 < k^{u,v) < . (18) 

V 

One also has k^(u,v)k_{u,v) = — . (19) 

Therefore, the equations (|16]) have a common root if and only if 

x^-" = K+(/?,5)K+(a,7). (20) 



Condition (|20D is an explicit constraint on the parameters of the ASEP. It 
defines the locus of points in the phase diagram accessible to n-dimensional 
representations. In Fig. 1 are shown the curves (hyperbolae branches) in 
the «:+(«, 7) — K+{P,6) plane where n = 1, ..,8 dimensional representations 
exist. It is worthwhile noticing that Eq. pO|) proves a conjecture made in 
which states that finite dimensional representations exist only in the region 

K+(/5,5)K+(a,7) > 1 . 

If the condition (|20| ) is fulfilled the common root of (|I6|) reads 
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and there exist vectors {W\ and \V) such that the boundary conditions 
are satisfied. These vectors can be computed as: 



{W\ 



k-1 



\v >-- 



V2 



m) 



with Vk 

and Wk 



n 



/3 [1 - X-™] [ax™ - K_(/3, 5)] 

n-2 _^ 

n ji 

7, 1 -L 



m=k 



[1 - [(ax™)-i - /t-(a,7)] 



for k = 1, 



(23) 
n . (24) 



With this definition, one has Vi = 1 and tu„ = 1. 

Hence, we have constructed an exphcit n-dimensional representation ( pUj , 
TT| , ^ and |24D of the quadratic algebra with the required boundary 
condrtions, provided that the rates a,/?, 7 and S satisfy the constraint (|20|). 
As discussed earher (section 2.1), we must verify that the matrix elemenfs~of 
length L (^, computed with this representation, do not identically vanish. 
In appendix B, using general properties of the algebra derived in appendix 
A, we prove that these matrix elements are non-trivial if the size L of the 
system is larger than the dimension n of the matrices. For L < n there are 
certain choices of the parameters for which, surprisingly enough, all matrix 
elements of length L vanish. We identify these choices of the parameters 
in appendix B. In any case, our representation can be used to investigate 
thermodynamic behaviour. 



3 Bulk density and correlations in the ther- 
modynamic limit 

In the thermodynamic limit the asymmetric exclusion process exhibits three 
different phases in which the current and correlation functions are given by 
different expressions. An exhaustive study has been carried out, mainly 
for the totally asymmetric model ^, |], O]. Much less is known about the 
Dartially asymmetric exclusion process: me phase diagram was obtained in 
M from an exact calculation of the current in the thermodynamic limit. This 
Diiase diagram, shown in Fig.l, agrees with the mean-field prediction. The 
Dulk density has been computed in Ref. H in mean- field approximation, 
and it was argued that the mean-field result is exact. We recall here the 
description of the phase diagram together with some known results. 
The phases can be described as follows: 
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Phase A (High Density): > K+(a,7) and 5) > 1. 

In the thermodynamic hmit, the current J a is 

J^ = (l-x) 

[l + /^+(A5)f 

and the mean-field prediction for the density in the bulk is: 

MF _ 

Pa 



1 + 5)- 

Phase B (Low Density): K+(a,7) > n+{j3,6) and K+(a,7) > 1. 

[1 + K+(a,7)] 

pr- ' 



1 + K+(a,7) ■ 

Phase C (Maximal Current): k,^{P,6) < 1 and «;+(«, 7) < 1. 



J, 



I — X 



c 



4 

Pd^ 



MF 1 



Phase A and B are separated by a coexistence line, defined by /t+(/3,5) = 
K+(a,7) > 1, through which the (mean-field) bulk density is discontinuous. 
This phenomenon has been investigated in the totally asymmetric case, and 
reveals the presence of a shock between a region of low density and a region 
of high density ||^, M. The same phenomenon happens in the partially asym- 
metric wilrbe confirmed by our calculation of correlation functions. 

Little is known about general correlation functions of the partially asym- 
metric exclusion process. In p|, some correlation functions were computed 
with one and two dimensional representations, i.e. on lines 1 and 2 in Fig.l. 
We shall now identify the generic structure of correlation functions for the 
case where n dimensional representations exist. Such an exact result can- 
not be obtained from a mean- field analysis |^. For example a two-point 
function is given by: 

= — — wmv) — - 

with C = D + E . (26) 
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The positions, on which a correlation function depends, enter above expres- 
sion as an exponent of the matrix C . For this reason we have to compute 
powers of the matrix C which is given by: 



C 



with \i 



D + E 



/ Ai 
1 A2 
1 A3 



2 + ax^-^ + 



ax 



k-l 



n—k 



\ 



(27) 



1 A„ J 



(28) 



We have used Eqs. (^0]) and (|2T|) to derive the last line. 

For L ^ 1, the dominating matrix elements of will be found in the 
invariant subspace of C that corresponds to the largest eigenvalue. We must 
identify the highest eigenvalue, and examine whether C can be diagonalized 
in the associated invariant spaceQ We notice that all the eigenvalues of C 
lie on the curve z t— > 2 + z + 1/ z , with z = a, ax, . . . , ax""^ . This implies 
that the largest eigenvalue of C can only be Ai, or A„ or both of them if it 
happens that Ai = An . We have: 



Ai - An = a + - - (ax"^ ^ + 

a \ ax 

= (l-a;"-i)(/t+(/3,5)-K+(a,7)) 



n-1 



[l-X 



1 



a — 



ax 



n-1 



(29) 



There are three different cases to consider, that naturally correspond to dif- 
ferent regions in the phase space: 



Bulk density and correlations in Phase A. This phase is defined by 
K+(/?,(5) > K+(a,7) and k+(/?, 5) > 1, therefore Ai > An- The leading 
eigenvalue of C is Ai and it is non degenerate. We write |Ai) (and (Ai|) for 
the corresponding right (and left) eigenvector. One should notice that (Ai| 
is equal to Hence we obtain that for large system sizes: 

Zl = ~ {w\x^){\i\v)\[ . (30) 

^We recall that the invariant space associated with the eigenvalue A is given by 
U^oKer(C-A)V 
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The bulk density is given by: 



PA = {n) = ^wicW) 

{W\\,)\\-'{X^\D\X,)Xt'{Xi\V) 
{W\Xi){Xi\V)Xf 

_ 1 + a 



(31) 
(32) 



Here we used that {Xi\D = (Ai|(l + a) and (Ai|Ai) = 1 to derive the last two 
lines. This result agrees with the mean-field prediction, as it was conjectured 
in i. 

Next we compute {W\C^-^DC''-^-'^DC^-''\V) . The result is a linear 
combination of powers oi the eigenvalues. Dividing it by the asymptotic 
expression (^) for Zl and using Eq. (|25| ) shows that two-point function is 
a superposition of terms (Aj/Aij™, where m is one of the distances j, k — 
j, L — k. Similarly, one can compute higher order correlation functions. For 
a large system, i.e. for L ^ n , they decay exponentially with correlation 
lengths 




" 2 + x-i«:+(/5,5) + a:"-«:+(a,7) j '-^'-'^ ' "^^^^ 

For n = 2 the correlation length computed in p|] is recovered. One should 
noticed that the correlation lengths clepend only on the functions k+(/?,(5) 
and K,+ {a, 7). 



Bulk density and correlations in Phase B. Phase B, defined by 7) > 
K+(P,6) and K+{a,-y) > 1 , is related to phase A by means of a symmetry: 
replacing Xj by 1 — r^+i^i as well as /3 by a, 5 by 7 and vice versa leaves the 
dynamics of the system invariant. We find that the bulk density is: 

PB = ^ (34) 



A, 



n 
1 



l + /€+(a,7) ^ ' 
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And the correlation lengths are given by: 

_ jln 2 + ^-'Mft^) + >^.K7) y 1.(36) 



Correlation functions on the coexistence line. Phases A and B coexist 
in a single system on the line k^(0,6) = K+(a,7) > 1. For the fully asym- 
metric process it was shown in g, g that the density profile depends linearly 
on the position, which indicates that the border between the two phases can 
be anywhere with the same probability. The same behaviour was observed 
in for the case where 2-dimensional representations exist, and two-point 
correlations functions were found to depend algebraically on the positions. 

Let us discuss the case where n dimensional representations exist. The 
eigenvalues of C are given by Eq. (EHI) which now reads 



Afc = 2 + [ x""-^ + x''-'' ] 6) (37) 

because fi;+(a,7) = . Obviously, these eigenvalues pairwise coincide: 

Afc = Xn+i~k for all k , and due to the line under the diagonal (see Eq. 
p7|) C is not diagonalizable. However, it can be transformed into a Jordan 
normal form with a two dimensional Jordan block for each eigenvalue. This 
implies algebraic behaviour of the correlation functions (see the discussion in 
Ref. |]13|)- We denote by (|Ai), \ \n)) a basis of the 2-dimensional invariant 
space associated with the highest eigenvalue A = Ai = A„ . In this basis the 
restriction of can be written as: 

(38) 

We choose |A„) = \n) and (Ai| = By decomposing on this normal basis, 
we find for large system sizes 

Zl = {W\C^\V) ~ LX^-^ , (39) 

because Wn = vi = 1. 

The m-point correlation function can now be computed as 
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with = xo < xi < 0:2 < . . . < < 1 • In the hmit of large system size, we 
project on the (|Ai), |A„)) plane and keep only highest order terms: 

( 1 L 2 L • • • m L ) 

k=l 

- E(^^-^^-i) ^ \: ^ — • (41) 
k=i ^ 

To obtain the last line we used that -D|A„) = (1 + ax"~^)|An) and (Ai|D = 
(Ai|(l + a) . Using furthermore the expressions (|3lD and ( p4D for the bulk 



densities in phases A and B we find: 

(r.,Lr.2L...r.„L) ^ Er=i(^fc - ^k-i)p'h-' p7''^' ■ (42) 

The formula ( ^21 ) clearly shows that on the coexistence line there is a shock 
between a low oensity phase ps and a high density phase pa- This shock can 
be anywhere between 1 and L with the same probability. 



4 Conclusion 

The finite dimensional Fock representations of the quadratic algebra, ob- 
tained here, allow us to derive exact results for the partially asymmetric 
exclusion process on some special curves of the phase diagram. We believe 
that the formulae obtained for the correlation functions are valid throughout 
phases A and B, but this is a conjecture: in order to prove it, one would have 
to use infinite dimensional representations, that would lead to complicated 
calculations. Maybe a kind of 'analytic continuation argument' could be used 
to obtain results vahd on any curve of the form 5)K-(.(a;, 7) = 

where n is not necessarily an integer any more, but we don't know how to 
do that. An interesting, and related, question is to understand the intuitive 
physics behind condition (^): what makes a system, that can be analyzed 
with finite n x n matrices, different from the others? 

The matrices that we have constructed are also useful in the representa- 
tion theory of the algebras associated with more general reaction-diffusion 
processes. In 0] the following system of algebraic relations is derived for such 
algebras: 

K1DE + K2ED = D + E, 

KsD^ = K4DE + K5ED, 

KqE'^ = KjDE + KsED, (43) 
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where the kj can be computed from the rates of the reaction-diffusion process 
[S]. The representations of (1431) are a subset of those of the deformed oscillator 
algebra. An analysis, similar to the one given in section 2.3, would give 
conditions on the rates for a finite dimensional representation of (^) to 
exist. 
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5 Appendix A 

In this appendix we show, without using any explicit representation, how to 
compute recursively expressions of the type 

{W\D'''E"'' ...D'"'E'^^\V) with rii + mi + . . . + Uk + ruk = L , 

that we shall call matrix elements of length L, and investigate conditions 
under which these quantities are different from zero. We need the following 
properties: 

Definiteness Property: Matrix elements of length L are either all strictly 
positive, or all strictly negative or all identically equal to zero. 

Indeed, matrix elements of length L are stationary solution of the master 
equation of the ASEP on a open chain with L sites [H. Therefore, according 
to the Perron- Frobenius theorem |^, they are unique up to a multiplicative 
constant and are all of the same sign. 



Reordering Property: The expression 

( . . . D"*^ E'^") - x" E^^i^i D^j=i , 

with L = ni + mi + . . . + Uk + mk and q = J2j<i<k ^jmi is equal to a linear 
combination of products of L — 1 operators D or E with strictly positive 
coefficients. 

One begins with the formula (obtained by induction on mi): 

jj^mi ^ a-rni^mi^ + (1 _ 3,) ^"^1 + y x'' E^ DE'^^-^-'' 

V ^0 
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that allows to transfer a matrix D from the left to the right of E"^^, and 
proves the reordering property for DE"^^ — x'^^E'^^D. One concludes then 
by induction on k and n^. 

The reordering property enables us to compute all matrix elements 
of length L if expressions of the type {W\E'-D^~'-\V) with I = 0, . . . , L and 
all matrix elements of length L — 1 are known. Hence one can relate matrix 
elements of length L to matrix elements of length L — 1. This can be achieved 
starting with the following system which is a simple consequence of (^): 

a{W\0-^^D^-^-^\V) - -f{W\D0D^-^-^\V) = {W\E^D^-^-^\V) . 
-S{W\E^D^-^-^E\V) + p{W\E^D^~^\V) = {W\E^D^-^-^\V) ' ^ ^ 

This system is rewritten according to the reordering property: 

a{W\E^+^D^-^-^\V) - x^-^{W\0D^~^\V) = Wl-i / 
-xL-i-^diWlE'+^D'^-'-^lV) + PiWlE^D'^-'lV) = ^ ^ 

where Wl-i and Wl-i are positive linear combinations of matrix element 
of length L — 1. We must distinguish two cases: 

(i) if aP — x^~^'y6 ^ 0, then all expressions of length L can be 
computed from expressions of length L — 1. Furthermore all matrix elements 
of length L are equal to if and only if all those of length L — 1 are equal to 
zero. 

(ii) If ajS — x^'^jS = 0, then all matrix elements of length less or 
equal to L — 1 identically vanish. Indeed, the l.h.s. of the equations in (^) 
are proportional, therefore one must have: 

= aW'L-i + x^-^-^6Wl-i (A3) 

Because of the definiteness property the positive linear combination on the 
r.h.s. can only be equal to 0, if all matrix elements of length L — 1 are equal 
to 0. For any n < L — 2, one has a(3 — x"''yS ^ 0, because otherwise this 
would imply a/3 = 7^ = which is impossible since all the rates are strictly 
positive. Using case (i), one concludes that all matrix elements of length less 
or equal to L — 1 are identically 0. 



a Wl-x 
-x^-^-^5 Wl-i 



In conclusion, we have proved the following properties of any repre- 
sentation of the algebra together with boundary conditions (Q): 
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Property A: If there is one matrix element of length L different from 0, 
then all matrix elements of length > L are non-zero. 

Property B: Suppose that there is one matrix element of length L different 
from 0, then for expressions of length < L, one has to consider two cases: 

(i) If for all non negative integer / a(3 — x^'jS ^ 0, then all the 
expressions {W\D"'^E"^^ . . . D'^'^ E^*'\V) are non-zero for any system size. 

(ii) If there exists one I such that a(3 — x^'^5 = 0, then all expressions 
of length < / identically vanish but all matrix elements of length > / are 
non- zero. 



6 Appendix B 

We first show, using Property A derived in Appendix A, that the represen- 
tations found in sections 2.2 and 2.3 provide non-trivial weights for systems 
of length larger than n or equal to n. Indeed the matrix elements (iy|D'|V^) 
for Z = 0, . . . , n — 1 can not all vanish, otherwise one would have: 

n 

= {W\D'\V) = J2{fx,yw,v, for / = 0,...,n-l, (Bl) 

1=1 

where /ij = 1 + ax*~^ is the ith eigenvalue of D. But the relations (IBl] ) 
can be interpreted as a system of n equations with n unknowns {wiVi). This 
system is a Van der Monde system, and as all the eigenvalues /ij are different 
from each other, the only solution is (wiVi) = for i = 1^ . , n. But this is 
not the case: as one can see from the explicit formulae and ^3), all the 
components of \V) and {W\ are different from 0. Consequently, the matrix 

elements (iy|-D'|V) cannot vanish for all /, i.e. there is an /' < n such that 

{W\D'- \V) 7^ 0. According to Property A, of appendix A, we conclude that 
the representation (|T0|) and (O), together with the boundary vectors given 
by (^2D and (^SF) always provides the stationary probabilities of the partially 
asymmetric exclusion process for systems of size bigger than n. 

Secondly, for systems of size strictly less than n, one has to check, 
according to Property B of appendix A, whether ajS — x^jS can be equal 
to for some non-negative integer k. Using (|T9|) and (|20|) , this condition is 
equivalent to: 

x"-!-'^ = K_(/3,(5)/t_(a,7). (B2) 
So, if (|B2|) is never satisfied, the n-dimensional representation found above 



allows to com pute the stationary weights for any system size. If there is a 



ko, such that (B2[) is true, then fco has to be less than n — 2 (because of |T3 



and the n-dimensional representation provides non-trivial weights only lor 
systems of size bigger or equal to ko + 1. 
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Caption to the figure 

Fig.l: Phase diagram for the ASEP with < x < 1 in terms of k+(q;, 7) and 
5) where v) = j^[—u + v + 1 ± ^{u — v — 1)"^ + iuv] . The phases 

are separated by the sohd hnes. The dashed hnes are the hnes where finite 
dimensional representations exist for a; = 3/4, and the numbers attached to 
these hnes are the dimensions of the corresponding representations. 
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